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Abstract 

We analyze the qualitative behaviors of scalar-tensor cosmologies with an 
arbitrary monotonic cj(<i>) function. In particular, we are interested on scalar- 
tensor theories distinguishable at early epochs from General Relativity (GR) 
but leading to predictions compatible with solar-system experiments. After 
extending the method developed by Lorentz-Petzold and Barrow, we establish 
the conditions required for convergence towards GR at i — ?■ oo. Then, we 
obtain all the asymptotic analytical solutions at early times which are possible 
in the framework of these theories. The subsequent qualitative evolution, 
from these asymptotic solutions until their later convergence towards GR, has 
been then analyzed by means of numerical computations. From this analysis, 
we have been able to establish a classification of the different qualitative 
behaviors of scalar-tensor cosmological models with an arbitrary monotonic 
cj(<i>) function. 

PACS number(s): 04.50.-Fh, 04.80.-Fz, 98.80.Cq, 98.80.Hw 
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I. INTRODUCTION 



The Einstein Equivalence Principle generates a whole class of theories (called metric 
theories) describing gravitational interaction. The distinction among these theories lies in 
the number and type (tensorial, vectorial, scalar) of additional gravitational fields, other 
than the metric tensor (/^j,, they contain. Erom a theoretical point of view, the most natural 
alternatives to General Relativity (GR) are scalar-tensor theories, which contain an addi- 
tional scalar field, (/>, the relative importance of which is determined by an arbitrary coupling 
function ^{(j)) [1-3]. In recent years, this class of metric theories has received a renewed in- 
terest specially in cosmology [4-11], because it provides a natural (non-fine-tuned) way to 
restore the original ideas of infiation while avoiding the cosmological difficulties coming from 
the vacuum-dominated exponential expansion obtained in GR. Scalar-tensor theories also 
arise in current theoretical attempts at deepening the connection between gravitation and 
the other interactions. Eor example, in modern revivals of the Kaluza-Klein theory and 
in supersymmetric theories with extra dimensions, one or several scalar fields arise in the 
compactification of these extra dimensions [12-18]. Eurthermore, scalar-tensor theories may 
also appear as a low-energy limit of superstring theories (see, e.g., Ref. [19]). 

Observational constraints upon scalar-tensor theories can be derived from their effects 
on solar-system experiments [20]. These quasi-stationary weak field tests imply that any 
alternative gravity theory must be, in that regime, very close to GR. However, these tests say 
nothing about how correct is a theory when gravitational forces are very strong. In general, 
the physical conditions in the early universe and the subsequent cosmological evolution in 
scalar-tensor theories could be very different from that obtained in the framework of GR. 
Since a wide class of these theories exhibits an attractive mechanism toward Einstein's 
theory [8,9,21-23], they could be compatible with present solar-system experiments in spite 
of their very different predictions at early times. The full consequences of an arbitrary 
function ^{(j)) will be only known when the cosmological solutions of the field equations are 
studied [3]. Analytical or numerical solutions for the scalar-tensor cosmological models are 
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well known in the framework of some particular theories proposed in the literature. This is 
the case of Brans-Dickers theory [24-29], Barker's constant-G theory [30,31], Bekenstein's 
variable rest mass theory [32-34], or Schmidt-Greiner-Heinz-MuUer's theory [35,36]. On the 
contrary, the properties of more general scalar-tensor theories still remain poorly known. 
Some important progress have been nevertheless made in the last years. For example, Burd 
and Coley [37] have used a dynamical system treatment to analyze the qualitative behavior 
of those models which result when a constant (Brans-Dicke) coupling function is perturbed 
by a slight dependence on the scalar field (e.g., ^{(j)) = cuq + </>~", with (/> — > oo and a > 0). 

A different approach has raised more recently from the work by Barrow [28]. In that work, 
he improved the method by Lorentz-Petzold [38] to solve, in an easy way, the cosmological 
equations of any scalar-tensor theory with a specified form of c<;((/)). By using this method, 
it has been possible to analyze the behavior at early times of some families of theories 
[28,39,40]. However, an exhaustive study of the properties of any possible scalar-tensor 
cosmology still is not available in the literature. The aim of this paper is to perform such a 
study in the hope of finding models with early behaviors qualitatively different from those 
found in the previous particular cases. To that end, we will restrict ourselves to theories 
with a vanishing cosmological constant term, and where ^{(j)) is a monotonic, but arbitrary, 
function of (j). 

The paper is arranged as follows. We begin outlining the Scalar- Tensor theories and 
we then show how to build up homogeneous and isotropic cosmological models in their 
framework (Sect. II). After establishing the conditions required for the convergence towards 
GR (Sect. Ill) we will present all the possible asymptotic solutions at early times which can 
be found in this class of theories (Sect. IV). The radiation-dominated evolution from these 
early behaviors until their later convergence towards GR will be obtained (Sect. V) for a 
wide class of scalar-tensor theories. Finally, conclusions and a summary of our results are 
given in Sec. VI. 
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II. SCALAR-TENSOR GRAVITY THEORIES 



A. Field equations: the Jordan and Einstein frames 

The most general action describing a massless scalar-tensor theory of gravitation is [1-3] 

s = -^ f{<pn- ^cp,^cp'^)^d*x + Sm (1) 

Idtt J (p 

where TZ is the curvature scalar of the metrics g^^^ g = det[g^^)^ (f) is the scalar field, and io{(f)) 
is an arbitrary coupling function determining the relative importance of the scalar field. The 
action (1) has been expressed in terms of the metric tensor (/^j,, called the "Jordan frame" 
[41,42], to which matter is universally coupled. In this frame, (/^j, is measured by using 
non-gravitational rods and clocks. 

The dynamics are sometimes more simply described by using the "Einstein frame" 
[41,42,23] 

S* = -fjT- I {n. - 2^,,^'^)y^^— + S*M (2) 
IdttG* j c 

where Newton's constant G* and the metrics g*^^ are now measured by using purely gravi- 
tational clocks. This frame is obtained from (1) by a conformal transformation 

A\^) = {G.cf>)-' (3) 

where A[(f) is an arbitrary function related to ^{(j)) by 

, 9 In A 

a= 3 + 2uj -'/^= ^— 4 
dip 

Since our measures are based on non-purely gravitational rods and clocks, observable 
quantities are those written in the Jordan frame. Comparison between theory and ob- 
servations must be then performed by using this physical frame, which will be that used 
throughout this paper. 

The variation of Eq. (1) with respect to g^^ and (j) leads to the field equations: 



- 9nu°(l^) (5a) 
(3 + 2cu)n</) = 8^r-cu'</),„</)'" (5b) 

which satisfy the usual conservation law 

T.f = (6) 

where T^" is the energy-momentum tensor, uj' denotes dio/dcj) and □(/> = g^^ (f^^fi-v 

B. Cosmological models 

In order to build up cosmological models, we consider a homogeneous and isotropic 
universe. The line-element has then a Robertson- Walker form: 



ds' = -c'dt' + R\t) 



1 - Kr^ 

and the energy-momentum tensor corresponds to that of a perfect fluid 



(7) 



T'^- = (p + p/c')u^u, + Pg^, (8) 

where K = 0, ±1, R{t) is the scale factor, p and P are the energy-mass density and pressure, 
respectively, and is the 4-velocity of the fluid. The field equations (5a) and (5b) then 
become 

Stt c^k Lli (j)^ acj) 

3</)'' -RqO^ 6 </)^ acj) 

(9) 

where a = R/Rq is the dimensionless scale factor and dots mean time derivatives. In 
addition, we have the conservation equation 

d{pa^) + {P/c'^)da^ = (10) 



Exact solutions of Eqs. (9) can be found, whatever the sign of 3 + 2a; is, by extending the 
method developed by Lorentz-Petzold [38] and Barrow [28]. During the vacuum [p = P = 0) 
and radiation-dominated epochs, the energy-mass density is given hy p = 3P/c^ oc a~'^. By 
introducing a conformal time, r, and the Lorentz-Petzold [38] variable, y, defined as 

adr = dt (11) 

y = (12) 

equations (9) then transform to 

y'2 = AVy + |</)'2a4(3 + 2u) (13) 
<^" + 2^</)' = -^^^g^ (14) 

where primes denote dj dr and V = Sirprad/'^ {Prad is the energy-mass density at the end of 
the radiation-dominated era). As usual, we have considered that k/a'^ can be neglected in 
the radiation-dominated field equations and, hence, that zero curvature models provide a 
good description of the early evolution of the universe. 
Integration of Eq. (14) gives 

</)'a2 = sign(3+2cu)2y3rA \ 3 + 2u \-^^^ (15) 

where A is a constant and where we have allowed for negative values of (3 -|- 2a;). 
By using Eq. (15), the integration of Eqs. (13) and (14) gives 

y = r[(r + Tof - sign(3 + 2u)A^] (16) 



anc 



l\3 + 2ur^ 



VSlnif) (3 + 2cu>0) 

2y3{arctan[(r + To) /A] - tt/2} (3 + 2cu < 0) 



:i7) 



where 



/ = ^^^ (18) 

To being an integration constant. 

Obviously, for a given form of c<;((/)), Eq. (17) can be integrated to obtain </>(r). Equations 
(12) and (16) give then a(r), and the relation r(t) can be finally obtained from (11). 

It must be noted that, from the relation a(r), it is also possible to obtain the speed- 
up factor, ^ = H/ H^^^ , where H = a/ a is the Hubble parameter, while i/™^ is that 
predicted by GR at the same temperature. To that end, we use the relation i/™^ oc 
together with Eq. (11), and we find 

^ oc a' (19) 

where ^ can be normalized so that >lasr— >-|-oo. 

Eor future purposes, it is convenient to write Eq. (17) in a more compact form. To that 
end, we define 



p = ^ln\y\ (20) 



and 



^=-i/|3 + 2.|'/^| (21) 

By using these quantities, Eq. (17) can be then written as (with K = — AF^^^ and 
C = sign(/f)^/2) 

-V^ln\Ke-P + (1 + /f2e-2p)i/2i (3 + 2cu > 0) 
07 = <; (22) 
-y3(arctan[ii^i|g^]-C) (3 + 2c. < 0) 



and, hence 



^ = -i|3 + 2c.r/^i^ (23a) 
dp 2 (f) dp 



VSKie^P + /f2)-i/2 (3 + 2cu > 0) 
-^3/^(62?' - /f2)-i/2 (3 + 2cu < 0) 



(23b) 



It is interesting to note that the function defined in Eq. (21) corresponds (except 
for a constant term) to the Einstein scalar field given by Eqs. (3) and (4). Similarly, Eq. 
(22) for 3 + 2a; > coincides with the radiation-dominated solution obtained by Damour 
& Nordtvedt [23] by using the Einstein frame. However, it is important to remark that 
we are always using the physical Jordan frame. We just use the quantities and p as an 
abbreviated notation allowing an easy comparison of some of our results and those obtained 
by using the Einstein frame. 

III. ASYMPTOTIC BEHAVIORS OF SCALAR-TENSOR THEORIES: 
CONVERGENCE TO GENERAL RELATIVITY 

We are only interested on viable scalar-tensor cosmological models, that is, those which 
are compatible with present solar-system observations but where the early evolution of the 
Universe is distinguishable from that obtained in GR. 

In order to ensure compatibility with solar-system experiments, we will later consider 
that any viable scalar-tensor gravity theory must converge towards GR (| |— > oo and 
|— > 0) as r — > -|-oo (equivalent to t — > oo). We show elsewhere [43] that primordial 
nucleosynthesis requires in fact that scalar-tensor theories which converge towards GR are 
"indistinguishable" from GR as soon as the end of the radiation-dominated era. 

A large bound on the present value of uj (hereafter loq) implies [3] that 

4 + 2cuo , , 

$0 = — ^ 1 24 

where $o is the present value of the dimensionless scalar field 

$ = Go<f). (25) 

Consequently, if we define x =| $ — 1 |, the asymptotic behavior of | 3 -|- 2a; | when r — > -|-oo 
(x — > 0) can be described as some power law 3X~^x~^ (with A > and e > 0). We have used 
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3A~^ for the proportionality constant to facilitate the comparison with the Eq. (53) below. 
It is known [39] that convergence towards GR ju^ |— > 0) requires e > 1/2. We will show 
now that compatibility with solar-system experiments also requires an upper bound on e. 
To that end, we will distinguish the case in which $ converges towards the GR value ($ = 1) 
as $ — > 1~, from the case in which $ — > l"*". 

If $ — > 1~ (x = 1 — $), the integral appearing in Eq. (17) and (21) can be then 
evaluated in a convenient closed form for a variety of values of e. When e/2 is a half-integer, 
we have 



4 lo 

2A 



1 



1 + 



(^-l)/2 9™-(2fc-l)/2' 

+ E 

k=i 



2k - 1 



(26) 



and when e/2 is an integer 



2A 



In 



1 — X 



(e-2)/2 

+ ^ V 

k=l , 



(27) 



In the same way, if $ — > l"*" (x = $ — 1), the integration of Eq. (21) for half-integral 
values of e/2 gives 



2A 



(-1)^ 



(^-l)/2 o -(2fc-l)/2] 

2arctan(y^)+ (-1)'+^ _ ^ 



f28) 



while, for integral values of e/2 



^3 e±2 
-^( — 1) 2 

2A ^ ' 



In 



X -\- 



(e-2)/2 



k+i ; 



k=i 



k 



(29) 



Erom Eq. (22), we see that — > as p — > -|-oo whatever the sign of (3 -|- 2a;) is. This 
condition can only be reconciled with x — > 0, as required by solar-system experiments, if e 
(integer or half-integer) is strictly smaller than 2 in the expressions (26)-(29). As a matter of 
fact, if e > 2, the solutions (26)-(29) are dominated by the terms with the form x~" ^ which 
diverge to infinity as x — > 0. It is important to notice that in the case analyzed by Barrow 
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[39], (e = 2, $ < 1 and 3 + 2a; > 0), the limit — > implies x — > 1/2. Barrow identified this 
solution with that obtained in GR because $ — >const ensures c<;($) — >const. However, we 
note that the constant value to which $ approaches is 1/2 and Eq. (24) then implies that 
the value to which uj approaches is not compatible with solar-system experiments. 

By generalizing the previous analytical results to non-integer and non-half-integer e values 
(we have tested the validity of such generalization by means of numerical computations), we 
then conclude that convergence towards GR as r — > -|-cx3, requires 

l/2<e<2 (30) 

IV. QUALITATIVE BEHAVIORS OF SCALAR-TENSOR THEORIES AT EARLY 

TIMES 

If $ becomes equal to unity [(f = 0) at a finite time, the theory is GR at all the 
radiation-dominated times, because in this case Eq. (22) implies K = and then (f is 
always vanishing. Moreover, in that case, Eq. (23a) implies that uj is infinite. However, an 
infinite uj is not equivalent to have GR, excepted if uj[^) is monotonic. Since we consider 
only monotonic uj[^) and we are only interested in models where the early evolution of the 
universe is distinguishable from that obtained in GR, we will consider that 3 -|- 2a; must 
be finite at r = 0. Its value can be however positive , vanishing or negative. We will now 
analyze separately these three possible cases by using the method described in Sect. II B. 

The criteria we have used to set the origin of r-times are the following: when models 
are singular, that is, when a(r) vanishes at some r, we take the constant Tq appearing in 
Eq. (16) so that we have the more recent singularity at r = 0. On the contrary, if the scale 
factor has some non-vanishing minimum, the origin of times will be taken so that there is 
only one minimum a- value for positive r values. 
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A. Possible values of <i> at r = 0. 



In order to represent adequately the early behavior of 3 + 2a; we need to know which are 
the possible values of $ at r = 0. We will now show that simple arguments considerably 
reduce the range of possible values of $(0)(= $(r = 0)). 

1. Case 3 + 2cj positive at t = 

If 3 + 2a; > at r = 0, $(0) must be +oo or 0. The other possibilities for $(0) are 
excluded. 

Since 3 + 2a;(r = 0) is finite, if $(0) diverges to — oo, then 3 + 2a;($) (as a function of $) 
has a nonvanishing horizontal asymptote. The early behavior of the coupling function could 
be then described by (3 + 2a;) ~| A; | at r — > 0, and Eqs. (17) would imply $ = —fV^. 
Consequently, = y$~^ < and the dimensionless scale factor would become imaginary. 

On the other hand, if $(0) has some nonvanishing finite value, we have a contradiction 
with the limit dip j dp — > \/3 as p — > — oo implied by Eq.(23b). As a matter of fact, in such 
a case we can choose Tq = A in Eq. (16) so that, from Eq.(12), the scale factor vanishes at 
r = 0. Erom Eq.(20), we then deduce p — > — oo as r — > 0. The resulting cosmological models 
are singular and $(p) has a nonvanishing horizontal asymptote at early times. Consequently, 
d<^ldp and, from Eq.(23a), d^p/dp 7^ ^3. 

2. Case 3 + 2a; vanishing at t = 

If 3 + 2cu(r = 0) ^ 0, then $(0) must diverge to ±00 (if 3 + 2cu 0" the case $(0) -00 
is excluded) or vanish. 

Erom the same arguments as before we can exclude any other possibility for $(0). In 
the case where $(0) diverges to —00 and 3 + 2a;(r = 0) — > 0~ (from negative values), Eq. 
(17) would lead to solutions with < . Nonvanishing finite values of $(0) would lead to 
values of d(f/dp at early times which are not compatible with that obtained from Eq. (23b). 
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3. Case 3 + 2cj negative at t = 



If 3 + 2ijj[T = 0) < 0, then $(0) must be positive and finite. 

The non-negative value of $(0) comes from the condition = y$~^ > where, from 
Eq. (16), y is always positive when 3 + 2a; < 0. On the other hand, $(0) must be finite 
because, otherwise, 3 + 2c<;($) would have a nonvanishing horizontal asymptote as $ — > oo. 
The early behavior of the coupling function could be then described by — (3 + 2a;) ~| A; | and 
Eq. (21) would imply (f oc ln($) — > oo as r — > 0, in contradiction with the finite value of (f 
implied by Eq. (22) when 3 + 2a; < 0. 
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B. Solutions at early times 

In order to analyze all possible early behaviors of scalar-tensor theories, the coupling 
function near r = will be expressed as a power law of (x =| $ — 1 |). Sometimes it 
will be necessary to add a constant k to get the assumed uj value at r = 0. This additive 
constant also allows an easier comparison with our analysis of section V. The negative power 
of X allows for larger convergence radii of our expressions. 

The analytical expressions of all the possible early behaviors are now listed. 



1. Case 3 + 2cj positive at t = 

• If3-|-2c<;>0as$— >-|-oo, the asymptotic behavior of the coupling function at early 
times can be described by 

3 + 2cu ~ (3/A2)(A; + x-1) (A > 0, A; > 0, x = $ - 1). (31) 

Integration of Eq. (17) then leads to 

27131 arctan ( ] - Vkln f V x-' + k -V k\ ^ ^yA^(o < A; < 1) (32) 

Wx-^ + k + Vk^j Wx-^ + k + Vkj ' ^ - ' 

(33) 



where / is given by (18). At early times [x — > +00), the dominant term in the left-hand 

side of Eqs. (32) is the second one, whatever the value of k is, and therefore 

1/2 

I / f\/Vk-l \ i + x/Vk 

a ~ 2JrGokT p p ^oc r(i+^/^'=)/2 ^ 

^ \1 + 2{2k - l)/VVfc + /2A/Vfc J 

\ r r r -2x1 \^ 

$ ~ 1 + -^{f + f~^' - 2) ^oc T-^l^^ oc ai+Vv^ (34) 

_ A/vT-l 

^ ^ OC r(^/^-i)/2 oc aVv^+i 



where we have again taken Tq = A. 
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We see that all models described by Eq. (34) are singular (a(0) = 0), but the ^ value 
at r = is: i) (^(0) — > -\-oo, if Xj^/k < 1; ii) (^(0) — > const., if Xj^/k = 1; or Hi) 
C(O) ^ 0, ifX/Vk > 1. 

This class of solutions is the most often studied in the literature [37,40]. It can be 
considered as a perturbed constant-cu Brans-Dicke theory. As in that particular theory, 
models have a power-law behavior at early times. However, since ^{(j)) has not been 
taken as a constant, the solar-system constraints on loq do not hold at early times. 
Consequently, this class of theories is specially well suited to be used in extended 
inflation models. Our analysis in terms of a, $ and ^ shows that, although all models 
have power law solutions at r — > 0, their speed-up factors can instead exhibit very 
different qualitative behaviors (see, e.g., the non-monotonic behavior of ^ in theories 
defined by Eq. (31) with X/Vk > 1). 

It is also important to note that the asymptotic solutions given by Eq. (34) represent 
a much wider family of theories than that directly defined by Eq. (31). Eor example, 
the case 3 -|- 2a; oc $^/($ — 1)^ studied by Mimoso and Wands [40] has a Taylor 
approximation at early times ($ — > -|-cx3) given by 3 -|- 2a; oc 1/2 -\- x"^, which is a 
particular case of Eq. (31). The solutions found by those authors are in fact contained 
in Eq. (34) in the limit of small r values. 

• If 3 -|- 2a; > as $(0) — > 0. We can describe this case by taking 

3 + 2cu ~ (3/A2)(A; + x-1) (A > 0, A; > -1, x = 1 - $). (35) 

Equation (17) then leads to 

VkTlln ( + ^ - a A ±1 \ _ 2y^arcsin(y^) = ln(/^)(-l < A; < 0) (36) 
\Vx~^ + k + \/k + 1 J 

\Vx-' + k + ./kTT ) \^x-' + k + -/k) ^ 

where, in the limit $ — > (x — > 1), the dominant term in the left-hand side is the first 
one, whatever the value of k is. Then 
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a ~ ^rG'o/(4(A; + l))r (/-(i+a/VI+T)(i + 2(2A; - 1)/^/^ + f ^/^)) 

$ ~ 4(A; + 1)/Wfc+T ^^^^ ^ ^T^^f^ (37) 

where we have taken Tq = A to get r = at the singularity. Arrows in the previous 
expressions denote the limit r — > 0. 

We see from these equations that the condition $ — > as r — > requires: «J < 
\j\/k + 1 < 1, or ii) 1 < \j\/k + 1 < 3. In the first case, models are singular and 
(^(0) — > +CX3 while, in the second case, models are nonsingular and (^(0) — > — oo. 

The particular case A; = of Eq. (35) has been previously considered by Barrow [39] 
and Mimoso and Wands [40], who found analytical solutions which are consistent with 
those given by Eq. (37). 

We have also performed some similar computations using other possible representations 
for uj{^) as, for example, 3 + 2a; = (3/A^)a;^^ The general solutions of Eq. (17) can 
then have a form which is different from that given by Eq. (37) but, in the limit r — > 0, 
all expressions become identical (a oc r^^"^)/^, $ oc r^, ^ oc r"^^"'"^)/^). 

2. Case 3 + 2cj vanishing at t = 
• If 3 + 2ijj[0) — > and $(0) — > +cx3, we can analyze this case by considering 

|3 + 2cu|~ (3/A2)x-^ (e>0,x = $-l). (38) 

The integral appearing in Eq. (17) has then the solutions given in Eqs. (28)- 
(29). However, since now x — > +(30, the dominant terms in these solutions are 
( — 1)(^"'""^)/^2 arctan(A/x) when e/2 is a half-integer, and ( — l)^^"'"^)/^ ln[a;/(a; -|- 1)] when 
e/2 is an integer in Eq. (28) and Eq. (29) respectively. It must be noted that only 
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the solutions of Eq. (17) in the case where e/2 is an half-integer satisfy the condition 
$ — > +(30 defining this class of asymptotic models. Consequently, the case where e/2 
is integer will be not included here. 

If 3 + 2u{0) 0+, Eqs. (17), (28) and (29) lead to 

a ~ ^TGo/2[{t + ToY - A^/' I cos(ln(/^/2^) ^ 
$ ~ cos-2[ln(/^/2^] ^ +00 (39) 
e ^ A^(e"/^ - l)/e"/2A = const. 

where / is given by Eq. (18), with Tq taken as Tq = A{e'^l^ + 1)/(1 - e'^l^) in order to 
have the more recent singular point at r = 0. 

Erom Eq. (39), we see that all these models are singular, with (^(0) — > const. A 
particular case of this class of behaviors is Barker's theory, defined by 3 + 2a; = 
(A2 = 3 and e = 1). 

On the other hand, if 3 + 2c<;(0) — > 0~, we obtain 

a = ^2rG'oA(l + [(r + ro)/A]2)i/2cos(^) ^oc r oc t^'^ 
$ = cos"^(5') ->(x r"^ oc a"^ (40) 
(A/y2)[l-cos(^/A)]i/2 

where g = A[arctan[(r + To)/A] — 7r/2] and we have taken Tq = Acot(7r/2A) (with 
A > 1/2) to satisfy the condition $(0) — > -\-oo (or cos((/) — > 0). This condition cannot 
be obtained for A < 1/2 and, hence, such a range of A values is excluded here. 

Models are again singular, with !^ — > const., where the constant is greater than unity 
for A > 1 and smaller than unity for A < 1. 

• If 3 + 2c<;(0) — > 0"*" and $(0) — > — oo, we can describe the asymptotic form of the 
coupling function by taking 

(3 + 2cu) ~ (3/A2)x-^ (e>0,x = l-$). (41) 
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The integral appearing in Eq. (17) has then the solutions given in Eqs. (26) and (27) 
but, since x — > oo, the dominant terms are now the logarithmic ones. For half-integral 
values of e/2, we obtain 



AS 2 _ a2 

e = i 



where we have taken Tq = in order to obtain r = at the singularity, and A = 1 in 
order to satisfy the condition $(0) — > — oo. 

We see that the expansion of the universe is equal (at any time) to that obtained in 
GR. This result can be understood by using Eqs. (42) and the field equations (9) to 
verify that the scale factor, a, then satisfies the standard Eriedman-Robertson- Walker 
equation, while the scalar-field satisfies an uncoupled differential equation. 

A similar behavior is found for integral values of e/2. 
• If 3 -|- 2ijj[0) — > and $(0) — > 0, we can represent this case by taking 

|3 + 2cu|~ (3/A2)(x-i - 1) (x = l-$). (43) 

Then, if 3 + 2cu ^ 0+, Eq. (17) implies 

-2arcsinyx = ln(/^) (44) 



and, then 



[{r + Tof - A'Y/' 



cos(ln/~^/2^ 
^oc cos-i(ln/-^/2) ^ oo 

$ = cos2(ln/-^/2) ^ (45) 

^ ^ oc - cos-2(ln /-^/2) ^ -oo 



17 



where Tq = A[e'^^^ + l)/(e'^/^ — 1) in order to have $(0) — > which defines this case. 

This choice of Tq imphes that a(0) — > +oo and models are then nonsingular, with a 
nonvanishing minimum value of a at Tmin — 2A/[e'^^^ ~ l)• 
If 3+2cu ^ 0", Eq. (17) gives 

— arcsin = A{arctan[(r + To)/A] — it/2} (46) 

and we get 



TG'o[(r + To)' + AY^' I cos(^) |-i^(x| cos(^) \-'^ +oo 
$ = cos^{g) (47) 
— > oc — COS~^((/) — > — oo 

where g = A[arctan(/) — 7r/2] and we have taken (for A > 1/2) Tq = Acot(7r/2A) in 
order to satisfy the condition $(0) — > 0. We can see from Eq. (47) that these models 
are again nonsingular. 

3. Case 3 + 2cj negative at t = 

• If 3 + 2cu(0) < and $(0) < 1, by taking 

-(3 + 2cu) ~ (3/A2)x-i (x = 1-$), (48) 

Eq. (17) leads to 



In 



1 - 



1 + 



2A arctan (r + tq) /A-tt (49) 



a = ^rG'o/4(r2 + A2)i/2l_i_£! ^ const 

4e^ 

$ = const < 1 (50) 

(l + ea)2 ^ ^ 

r (l + e^)r 2 
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where g = A[arctan[(r + Tq) j A) — 7r/2] with Tq = 0. 

We see that a and $ never vanish and models are then nonsingular. 

• If 3 + 2cu(0) < and $(0) > 1, we can take -(3 + 2u) ~ (3/A2)a;-^ with x =| $ - 1 1 = 
$ — 1. Solutions are then those given in Eqs. (40) but, now, with A < 1/2 in order to 
have a finite $(0). In this case, cos((/) never vanishes and models are then nonsingular. 

C. Existence of an initial singularity 

The existence of an initial singularity can be characterized in terms of the sign of 3+2c<;(0). 

When 3 + 2c<;(0) is negative or vanishing (— > 0~), we know that $(0) cannot diverge to 
infinity. Since the variable y never vanishes when 3 + 2c<;(0) < 0, Eq. (12) implies that all 
models are then nonsingular. The analysis of section IV B 3 shows that, in this case, the 
speed-up factor !^ is negative at r = 0. The negative value of !^ means a contraction phase 
of the Universe. 

On the contrary, if 3-|-2c<;(0) is positive or vanishing (— > 0"*"), we can choose Tq in Eq. (16) 
so that y oc r — > at early times. Erom Eq. (12), we then find that only those models where 
$(0)r" — > (a > 1) can be nonsingular. Erom Eqs. (45) and (37) we see that the condition 
a > 1 is always satisfied when 3 + 2c<;(0) — > 0"*" and $(0) — > 0. However, if 3 + 2c<;(0) — > 
const, that condition requires 1 < \j\/k + 1 < 3. As before, sections IV B 1 and IV B 2 show 
that the speed-up factor ^ initially diverges to — oo. 

In all the other cases, models are singular and sections IV B 1 and IV B 2 show that the 
initial value of the speed-up factor ^ can be 0, a positive constant, or -|-oo. 

Table 1 summarizes all the asymptotic solutions obtained in this Section as well as some 
of their main properties. Singular and non-singular models are denoted in the last column 
of this table by S and NS, respectively. 
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V. RADIATION-DOMINATED EVOLUTION IN SCALAR-TENSOR 

COSMOLOGICAL MODELS 



The evolution from the above asymptotic behaviors until their later convergence towards 
GR, can only be known by specifying a particular c<;($), assumed to be valid during all 
that time interval. In order to get some insight into the possible evolution of scalar-tensor 
cosmological models, we will write c<;($) as the following Taylor series. 

1 



I I 8 = 



^a, (51) 



the first order of which i 



IS 



ai|$-l| (52) 



where, in order to allow for convergence towards GR, we have taken Oq = and 1/2 < e < 2 
(see Sect. III). 

Eq. (52) is not just a generalization of that proposed by Garcia-Bellido and Quiros [22]. 
It gives an exact representation for most of the particular scalar-tensor theories proposed 
in the literature [44] and, what is more, it is also a first-order approximation to any other 
theory provided that $ is not very different from unity. Solar-system experiments imply 
that this last condition is in fact satisfied from some time up to the present. For t < t*, 
Eq. (52) then represents a particular class of theories which allows nevertheless to reproduce 
all the asymptotic behaviors described in the previous section. 

A monotonic c<;($) function, as that given by Eq. (52), implies that: i) $ is never equal 
to unity except for p — > -|-cx3, and ii) $ must have a monotonic time evolution. As a matter 
of fact, if $ becomes equal to unity at a finite time, the theory is GR (i.e. (f = 0) at that 
time. Since, in that case, Eq. (22) implies K = 0, such a theory then reduces to GR at all 
the radiation-dominated times. On the other hand, if $ is not a monotonic function of time, 
Eq. (23) would imply an infinite value of uj when / dp vanishes and, in turn, this would 
imply $ = 1. The theory then becomes identical to GR. Einally, since $ cannot diverge to 



20 



infinity except at r = (because we then have = y$ ^ — > and this defines the time 
r = 0), $ is confined to values $ > 1, if the present value of $ (hereafter $o) l"*" or to 
$ < 1, if $0 ^ 

In order to depict all possible cosmological models, we have then performed a numerical 
integration of the field equations (9) by using Eq. (52) and the temperature T as variable. 
These numerical solutions can be easily interpreted by considering the asymptotic behaviors 
found in previous sections. In order to facilitate the discussion, it is convenient to rewrite 
the function 3 + 2a; appearing in Eq. (17) as 

\3 + 2uj\= {3/X'^){x-' + k) (53) 

where A = ^3 |ai|^ /2, k = (AV3)(3 + 26), and 

1 - $ (if $ < 1) 

(54) 

$ - 1 (if $ > 1) 

where the condition $o < 1 is obtained for ai < and 3 + 2c<;o > 0, or ai > and 3 + 2c<;o < 0, 
while $0 > 1 is obtained for ai > and 3 + 2ijJo > 0, or ai < and 3 + 2ijJo < 0. 

We will now discuss the different cosmological models in terms of the k and A parameters. 
The qualitative behaviors of these models does not depend on e provided that 1/2 < e < 2. 



A. Models with cjq > 

1. Case <i>o —7- 1~ and k > —1 

In this case, Eq. (53) implies that 3 + 2a; is strictly positive at any time and that $ 
remains always smaller than unity. According to that quoted in section IV A, $ vanishes at 
r = 0. The early evolution of the universe then corresponds to the asymptotic solutions of 
section IV B 1 (Eq. (37)) and depends on the X/\/k + 1 value. 

• If < Xj^/k + 1 < 1, cosmological models are singular. ^ decreases monotonically 
from +00 at r = to its present value — > l"*" (Eigure 1). 
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• If 1 < \j\/k + 1 < 3, cosmological models are nonsingular. The scale factor has a 
minimum value at some time T^in^ where = 0, and the temperature of the universe 
has then a maximum value {T[Tmin) (Figure 2). 

If A; = and e = 1, the neglected terms in Eq. (36) are exactly vanishing. Hence, the 
analytical solutions given by Eqs. (37) are exact at any time. We have therefore found 
nonsingular models which, unlike those obtained by Barrow [39], converge towards GR 
as r — > +00 and, hence, they are compatible with all solar-system experiments. 

• If A; = and A = 1, solutions are those given by Eq. (42) and the scalar field is then 
uncoupled from the expansion of the Universe is then equal to unity at any time). 

2. Case $0 1"*" and k > 

In this case, Eq. (53) again implies that 3 + 2a; is strictly positive at any time and 
$(0) — > +00. The early evolution of these models then corresponds to the asymptotic 
solutions given by Eq. (34). 

Models are singular, with $ monotonically decreasing from $(0) — > +00. The behavior 
of depends instead on the value of Xj^/k: 

• If \j\fk < 1, 1^ decreases monotonically at any time until reaching unity (Eigure 3). 

• If \j\fk = 1, (^(0) — > const. Numerical integration shows in this case that the finite 
value of (^(0) is always greater than 1 (Eigure 4). 

• If Xj^/k > 1, then initially increases from (^(0) — > and later reaches a maximum 
value ^max which can be arbitrarily high. The ^max value increases with A;, while it 
slightly decreases with A. The radiation-dominated behavior of ^(t) is then highly 
non-monotonic, with an initial phase where ^ is slower than in GR and a later phase 
with a faster expansion of the universe which converges towards the GR value as 
T — > -|-oo (Eigure 5). No similar behaviors have been found for any other choice of 
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parameters, even for negative loq values. We have not either found models where ^ has 
some minimum value at r > 0. 

3. Case $0 1~ and k = -1 

3-\-2ijj is always positive and vanishes at r — > 0, with $(0) — > 0. The radiation-dominated 
evolution is then described by Eq. (45), which are exact solutions for e = 1, and models are 
therefore nonsingular. The qualitative behavior of these models is then similar to that found 
in section V A 1 when 1 < Xj^/k + 1 < 3. However, in this case, the speed-up has a smooth 
maximum value greater than unity after of which it falls quickly to zero (Figure 6). Note 
that solutions given by Eq. (45) provide another class of nonsingular models compatible 
with solar-system experiments. 

4. Case <i>o —7- l"*" and k = 

3 -|- 2a; is always positive and vanishes at r — > 0, with $ — > -|-oo. Solutions at any time 
are those given by Eqs. (28) and (29) while, at early times, they are approximated by Eqs. 
(39). Their early behavior is then similar to that discussed in section VA2 [Xj^/k = 1) 
(Eigure 4). These models are then singular, with $(0) — > -|-cx3 and ^ — > const > 1. 

5. Cases <i>o — ?■ 1~ and A; < —1 or <i>o —t- l"*" and k < 

These values of the parameters in Eq. (53) imply that 3 -|- 2a; vanishes at some finite 
temperature, where $ = Eq. (23) then implies that / dp — > — 00 as $ — > and, from 
Eqs. (12) and (19), we have also that ^ — > -|-oo in this limit. This class of theories is very 
unsatisfactory from a theoretical point of view. As a matter of fact, if the sign of 3 -|- 2a; 
changes at a finite temperature, it is not possible to set the integration constants in Eqs. 
(22) and (23) to avoid a discontinuity in (f or d(f/dp. Numerical integration produces in fact 
overfiows at this point (because of the infinite values of a', $' and (3 -|- 2a;)~^). We have not 
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tried to handle this hmit because we consider that such a change in the sign of 3 + 2a; is not 
physicaL Numerical integration stops then at such a temperature (see Figure 7). 



B. Models with cjq < 

1. Case $0 1~ and A; < 1 or $0 ^ 1"*" and k < 

In the case $0 1~ and k < 1 ($0 l"*" and k < 0), $(0) has some nonvanishing value 
smaller (greater) than unity and 3 + 2a; grows (decreases) converging towards (3/A^)(A; — | 
1 — $(0) |~^) at r — > 0. Numerical integrations show however that $(0) is very small (large) 
and, therefore, 3 + 2a;(0) has some nonvanishing negative value. The early behavior is then 
similar to that discussed in section IV B 3 (Eq. [50]) for $ < 1 ($ > 1). Consequently, 
models are nonsingular and their qualitative behavior is similar to that shown in Figure 2. 

2. Case $0 1~ and k = 1 

In this case, $(0) — > 0, 3 + 2a; is negative at any time, and it approaches 0~ at r = 0. 
The asymptotic behavior of this model is given by Fq. (47). As before, it is nonsingular and 
the evolution of ^ is that of Figure 2. 

3. Case $0 1"*" and k = 

In this case, $(0) — > +cx3, 3 + 2a; is negative at any time, and it approaches zero as 
r — > 0. The early evolution of ^ is given by Fq. (40) and depends on the A parameter: 

• If A > 1, then ^(0) > 1, and the qualitative behavior is similar to that of section V A2 
for Xj^/k = 1 (see Figure 4). 

• If A = 1, then ^(0) = 1. Numerical integration shows that ^ = 1 at any time and, there- 
fore, the scalar-field is uncoupled from the expansion of the universe. The behavior of 
these models is then equal to that of GR (see section V A 1). 
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• If A < 1, then ^(0) < 1, and the quahtative behavior is that shown in Figure 8. These 
are the only models we have found with a non- vanishing value of ^(0) smaller than 1. 

4- Cases <i>o — ?■ 1~ and k > 1 or <i>o —t- l"*" and k > 

These models imply that 3 + 2a; vanishes at some finite temperature. Their behavior has 
been described in section VA5 (Figure 7). 

VI. CONCLUSIONS 

We have analytically and numerically studied all the possible behaviors of scalar-tensor 
cosmologies with an arbitrary monotonic c<;($) function. Our main results are: 

1. The convergence towards GR during radiation-dominated era, requires an asymptotic 
behavior of the coupling function as | 3 -|- 2a; |oc| 1 — $ |~^, with 1/2 < e < 2. On the 
other hand, Damour & Nordtvedt [23] have found, by using the Finstein frame, a natural 
attractor mechanism toward GR during the matter-dominated era. Although the nature of 
these two conditions are totally different, it is straightforward to show that our condition is 
mathematically equivalent to that of Damour & Nordtvedt. Consequently, any model which 
does not satisfy our condition on the e parameter, will not either converge toward GR during 
the matter-radiation era. 

2. We have shown that the value of the $ field at r = and the existence of an initial 
singularity are constrained by the sign of 3 -|- 2a; at that time. In particular, if 3 -|- 2a; > 0, 
$(0) cannot have a nonvanishing finite value. On the contrary, if 3 -|- 2a; < 0, $ must be 
finite and all the resulting models are non-singular. 

3. The analytical asymptotic solutions of all the possible scalar-tensor cosmological mod- 
els converging towards GR have been discussed in terms of the sign of 3-|-2a; at r = 0. These 
solutions describe all the early behaviors for any monotonic form of the coupling function 
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Although some of these solutions are known in the literature^ (e-g-, [28,39,40]), we 
have been able to find several other behaviors which are qualitatively very different from the 
previous ones. We remark in particular that several of these solutions imply non-singular 
models which, unlike that found by Barrow [39], are compatible with solar-system experi- 
ments. 

4. From a specified form of the coupling function c<;($), which allows nevertheless to 
reproduce all the above asymptotic behaviors, we have numerically determined the evolution 
of the universe at any time. Our analysis reveals the existence of scalar-tensor cosmological 
models with qualitative behaviors very different from those previously found in the particular 
cases studied in the literature. We summarize in Figures 9 all these behaviors in terms of the 
A and k parameters defined in Fq. (53). It is important to note that, for the first time, we 
have found a class of models with an explicit non-monotonic time evolution of the speed-up 
factor. 

Finally, additional constraints on all these behaviors can be obtained from the light 
elements primordial production. This question is examined in a forthcoming paper. 
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FIGURES 

FIG. 1. Singular models with ^^(O) = and ^ decreasing monotonically from infinity (A^ = 3/4, 
k = 0, solid line corresponds to loq = 10^^ and dashed line to loq = 10^'*). Note that the largest loq 
values the most shifted to the left are the curves. However, the form of curves does not depend on 

FIG. 2. Non-singular models with ^ increasing from (A^ = 3/2, A; = 0, solid line corresponds 
to LJo = lO^-^ and dashed line to lvq = 10^^) 

FIG. 3. Singular models with ^^(O) — ?■ +00 and ^ decreasing monotonically A^ = 3/2, k = 26.5, 
solid line corresponds to loq = 10^° and dashed line to loq = 10^^) 

FIG. 4. Singular models with ^^(O) — ?■ +00 and ^ decreasing from a finite constant value, 
^(0) > 1 (A^ = 3/2, A; = 3/2, solid line corresponds to ujq = 10^® and dashed line to ujq = 10^°). 

FIG. 5. Singular models with ^^(O) — ?■ +00 and a highly non-monotonic ^ (A^ = 3/2, k = 1/2, 
solid line corresponds to ujq = 10^® and dashed line to ujq = 10^°) 

FIG. 6. Nonsingular models with ^^(O) = and ^ has a smooth maximum value greater than 
unity (A^ = 3/4, k = —1, solid line corresponds to ujq = 10^"^ and dashed line to ujq = 10^'*) 

FIG. 7. Models where the sign of 3-|-2cj changes at a finite temperature, T*. On this figure we 
plot a model with ^(r*) > 1 and ^ decreasing monotonically from infinity (A^ = 3/2, k = —2.5, 
solid line corresponds to ujq = 10^^ and dashed line to ujq = 10^"^). The opposite behavior is also 
possible with ^(r*) < 1 and ^ increasing monotonically from —00. 

FIG. 8. Singular models with ^^(O) — ?■ +00 and ^ increasing monotonically from a non-vanishing 
value of ^(0) < 1 (A^ = 3/4, A; = 0, solid line corresponds to loq = —10^^ and dashed line to 

LOO = -1023) 
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FIG. 9. Space of theories for a) cjq > 0, $0 < 1 (solid and dashed lines correspond to A^/(A;+1) = 
1 and 9, respectively), b) cjq > 0, $0 > 1 (the line corresponds to A^A; = 1), c) ujq < 0, $0 < 1, and 
d) LJo < 0, $0 > 1- 
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Table 1 



Summary of asymptotic solutions 



3 + 2cu(0) $(0) 


3 + 2cu ~ 


Solution 






> 


+ oo 


(3/A2)(A; + x-i) (k > 0) 


Eq. (34) 


+ 00 {X/^/k < 1) 


S 










const. [Xj^/k = 1) 


S 










( x/Vk > 1) 
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> 





(3/A2)(A; + x-i) (A; > - 


1) Eq. (37) 


+00 {X/Vk + 1 < 1) 
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( 1 < X/Vk + 1 < 3) NS 





+ 00 


(3/A2)x-^ (e > 0) 


Eq. (39) (0+ ) const. 


S 








Eq. (40) (0-) 


const. (A > 1/2) 
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-oo 


(3/A2)x-^ (e > 0) 


Eq. (42) 


1 (GR-like) 
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(3/A2)(x-i-l) 
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-oo 
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< 1 


-(3/A2)(x-i-l) 
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-oo 


NS 
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> 1 


-(3/A2)(x-i-l) 
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-oo (A < 1/2) 
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Fig.7 
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Fig.8 
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